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Superconductivity and Spin gap in the zigzag chain t-J model simulating a CuO 

double chain in P^Ba^urOis-^ 
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Using the numerical diagonalization method, we examine the one- dimensional t\-ti-J\-J% 
model (zigzag chain t-J model) which represents an effective model for metallic CuO double 
chain in the superconductor Pr2Ba4Cu70is_5. Based on the Tomonaga-Luttinger liquid theory, 
we calculate the Luttinger-liquid parameter K p as a function of electron density n. It is found 
that superconductivity is realized in parameter region corresponding to the experimental result. 
We show phase diagram of spin gap on the £2/^1 \-n plane by analyzing the expectation value of 
twist-operator Z a in the spin sector. The spin gap appears in the region with large £2/|£i |, where 
the phase boundary at half-filling is consistent with that of the known frustrated quantum spin 
system. The analysis also suggests that the estimated value of the spin gap reaches ~ 100K in 
the realistic parameter region of Pr2Ba4Cu70i5_^. 

KEYWORDS: zigzag chain t-J model, spin gap, numerical diagonalization, superconductivity, Tomonaga-Luttinger 
liquid theory 



Recently, Matsukawa et al. have discovered a new 
superconductor Pr2Ba4Cu70i5_a(Pr247) in which CuO 
double chains are considered to derive the superconduc- 
tivity at T c ~ 20K. 1 ' 2 Since electronic conduction in 
Cu02 plane of Pr247 is suppressed due to the so-called 
Fehrenbacher-Rice state, 3 the double chains are expected 
to play a crucial role for metallic state of the material. In 
fact, anisotropy in the resistivity of a single crystal shows 
the one-dimensional(lD) conductivity based on the CuO 
double chains and the NQR experiment also indicates 
that the superconductivity is realized in the CuO dou- 
ble chains. 4-7 These experiments stimulate our interest 
in the theoretical aspect for the electronic state and the 
superconductivity of the double chain system. 

Many theoretical works have been performed on the 
electronic state of double chain systems such as Ladder 
models and zigzag chain models. 8-23 Generally speak- 
ing, the electronic state of these two-band models is 
characterized by existence of four Fermi points, namely, 
±kp 1 and ±kp 2 on the Fermi surface. In the weak cou- 
pling regime, bosonaization method reveals that the low- 
energy excitations of the double chain are given by a 
single gapless charge mode with a gapped spin mode (la- 
beled as clsO), when the ratio of the two Fermi veloci- 
ties \vf 1 /vf 2 \ is smaller than a critical value ~ 8.6. 8-10 ' 24 
The correlation functions of the superconductivity(SC) 
and that of the charge density wave (CDW) decay as 

~ r and ~ cos[2(fc Fa — k Fl )r]r~ 2K " , respectively, 
while that of the spin density wave (SDW) decays ex- 
ponentially. Here, K p is the Luttinger liquid parameter 
and determines the critical exponents of various types of 
correlation functions in the model which is isotropic in 
spin space. 11-13 In the clsO region, the SC correlation is 
dominant for K p > 1/2, while, the CDW correlation is 
dominant for K p < 1/2. 

In the strong coupling regime, the double chain sys- 
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terns have been studied by using numerical methods. 14-21 
At half-filling, the system can be described by a Heisen- 
berg model whose ground state is a spin liquid insulator 
with a finite spin gap. 15, 17,21 Away from half-filling, the 
system becomes a metallic state which maintains a spin 
gap. 16,21 This behavior is explained by the existence 
of electron pairs produced by the dominant fluctuations 
of the Akp charge density wave or the interchain paring 
fluctuations. 

Among them, there are few works which consider the 
model just corresponding to the Pr247 except our pre- 
vious work 22 and the very recent work using the fluc- 
tuation exchange (FLEX) approximation. 23 In our pre- 
vious paper, we have investigated the superconductivity 
in the d-p double chain model, simulating a CuO dou- 
ble chain of Pr247 where the tight-binding parameters 
are determined so as to fit the band structure of the 
local density approximation(LDA). On the basis of the 
Tomonaga-Luttinger liquid theory, we have obtained K p 
as a function of the electron density n. The doping de- 
pendence of K p is in good agreement with that of T c 
in Pr247 2 when we assume that T c is a monotonically 
increasing function of K p at K p > 1/2. However the 
Hartree-Fock(HF) approximation has been used in this 
work and the analysis is limited in the case of the weak 
coupling region as well as the FLEX approximation. 

Since the strong correlation effect may play an impor- 
tant role in the electronic state and the superconductiv- 
ity of Pr247, a nonperturbative and reliable approach is 
required. In this work, we employ the numerical diago- 
nalization method for the double chain t-J model whose 
parameters are selected to cover the realistic band struc- 
ture of the CuO double chains. We calculate the Lut- 
tinger liquid parameters K p and address the behavior of 
the spin gap in the strong coupling regime beyond the 
previous works. 

We consider the following Hamiltonian for the one- 
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dimensional t\-t 2 -J\-J 2 model(zigzag chain t-J model); 
H = ti^(c\ tCJ c l+l . a + h.c.) + t 2 ^2(cl a c i+ 2,a + h.c.) 



t 2 ,J2 



* *- 



+ Ji y^(S t ■ S t+ i - -njn t+ i) 



where cj a stands for the creation operator of an electron 

with spin a at site i and nj i<7 = cj^c^. Here, t\ is the 
hopping energy between the nearest-neighbor sites and t 2 
is that between the next nearest-neighbor sites as shown 
in Fig. 1(a). The interaction parameters J\ and J 2 stand 
for the exchange coupling between the nearest-neighbor 
sites and between the next nearest-neighbor sites, respec- 
tively. We take account of the infinite on-site repulsion 
by removing states with doubly occupied sites from the 
Hilbcrt space. 

To determine the hopping energies t\ and t 2 , we com- 
pare the noninteracting t\-t 2 band with the d-band of 
the CuO double chain obtained by the LDA band struc- 
ture in YaBa 2 Cu 4 8 (YBCO). 25 Here YBCO includes 
the CuO double chains with the same lattice structure 
as those in Pr247. As shown in Fig. 1(b), both bands 
arc in good agreement with each other, when we select 
h = -O.lcV and t 2 = -0.45cV. 26 

As for the exchange coupling energies, Ji is consid- 
ered to originate in the exchange interaction between 
electrons in the nearest neighbor d-sites. The value of 
Ji is given by the 2nd order perturbation with respect 
to the hopping tdd, i-e. J\ — At 2 dd /Ud, where Ud is the 
on-site Coulomb interaction between d electrons. When 
we assume tdd = 0.12eV 26 and Ud = 6eV, Ji is estimated 
as O.OleV. On the other hand, J 2 is considered to the 
superexchange interaction acting between the next near- 
est neighbor Cu-sites connecting through O-site. Our 
previous study for d-p single chain and Ladder models 27 
indicates that the exchange interaction J/t p d is <~ 0.14 
for d-p single chain, and about 0.09 for d-p Ladder at 
A/tpd <~ 2.6 and Ud/t p d = 8. Here, we note that the 4th 
order perturbation 28 with respect to t p d overestimates 
the value of Jjt v d in the case of A/t p d ~ 4. 27 By ref- 
erence from the above, we regard the adequate value of 
J 2 as <~ 0.15cV. These values are close to that of corre- 
sponding parameters obtained in Cu0 2 plain system. 29 

We numerically diagonalizc the Hamiltonian up to 
24 sites using the standard Lanczos algorithm and cal- 
culate the ground state energy Eq. We use the pe- 
riodic boundary condition for N e = Am + 2 and the 
antiperiodic boundary condition for iV e = 4m, where 
N e is the total number of electrons and m is an inte- 
ger. The filling n is defined by n = N e /N, where N 
is the total number of sites. The critical exponent K p 
is related to the charge susceptibility Xc and the Drude 
weight D by K p = ^XcD) 1 ' 2 , with D = 
where Eo((fi) is the total energy of the ground state as a 
function of magnetic flux N(f). 13 Here, the flux is im- 
posed by introducing the following gauge transforma- 
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Fig. 1. (a) Schematic diagram of t\ — (2 — Ji — J2 model on the 
zigzag chain, (b) Energy dispersion relation for the noninteract- 
ing t\ — t2 model on the zigzag chain. Solid lines are the tight- 
binding result with t\ = — O.leV and t2 = — 0.45cV. Closed 
circles are the LDA result for the c(-band of the CuO double 
chain. 
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Fig. 2. K p as a function of n for ti = —O.lcV, (2 = — 0.45cV, 
J\ = O.OleV and J2 = 0.15eV. The dashed line is a guide 
for eyes. Inset shows the chemical potential fi as a function of 
n. The solid line represents a fitting line using a second order 
polynomial by the least square method. 



tion: c mcr f — > e lm< At a for an arbitrary site m. When 
the charge gap vanishes in the thermodynamic limit, 
the uniform charge susceptibility \c is obtained from 
Xc = §77, where the chemical potential n(N e ,N) is dc- 
fined by M (7V e ,A) = W±MQ^feiA) , Using the 



above Xc and D, we calculate the K p from the ground 
state energy Eq of the finite size system. 

Fig. 2 shows the Luttinger parameter K p as a function 
of the electron density n for t\ = —O.lcV, t 2 = — 0.45eV, 
J i = O.OleV and J 2 = 0.15eV. Inset shows the chem- 
ical potential /x as a function of n, where data of \i is 
fitted to a second-order polynomial as a function of n 
by the least square method and the value of Xc is esti- 
mated from differential coefficient of the polynomial. K p 
increases with increasing n and then have a maximum at 
an optimal electron density at n ~ 0.6. In the region of 
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Fig. 3. (a) K p as a function of J2 for n = 2/3(12clcctrons/18sites). 
(b) The energy difference Eo(<j>) — Eo(0) as a function of an 
external flux eft for n = 2/3(12clcctrons/18sitcs) at J2=0, 0.15 
and 0.30eV with ti = -O.lcV, t 2 = -0.45cV and Ji = O.OleV. 





ar' 


j 1 1 1 1 1 1 1 1 
t 1= 1 J 1= 0.4 


n=0.5/'' 




J 2 =0.4t 2 2 












=06. 

n=0.667 






n=0.75 # " 


..-••-"*" 










n=1 . 


22 20 


12" 


10 N=8- 

1 1 1 


) 0.005 




0.01 0.015 



1/N 2 

Fig. 4. Size dependence of the critical point *2 C for n = 
1/2,3/5,3/4,4/5, and 1 determined by the condition (Z a ) = 0. 
The open square stands the well known result of the J1-J2 
Hcisenberg model in the limit N — » 00. 



the other correlations (SC phase). The overall behavior 
of K p is consistent with our previous work obtained by 
the HF approximation 22 except near the half-filling. 

Figure 3(a) shows the value of K p as a function of 
J 2 for the 12clcctrons/18sites system. As J increases, 
K p increases and it becomes larger than 0.5 even if 
J 2 = 0. To confirm the superconductivity, we calcu- 
late the ground state energy E (<f)), as a function of an 
external flux <j>. As shown in Fig. 2(b), anomalous flux 
quantization clearly occurs at J > 0.05. It suggests that 
the SC phase appears at K p > 0.5. 

Next, we consider the phase diagram of the spin gap 
A CT . Generally speaking, it is not easy for numerical 
methods to estimate A CT preciously in the case of the 
energy scale being small. Especially, it is very difficult 
to determine the phase boundary of the spin gap which 
is defined by = 0. To overcome this difficulty, we 
introduce twist-operator Z a , given as 

2ni N 

Za = exp[— ^i(n jT - n n )]. (I) 

When the expectation value (Z a ) > 0(< 0), the system 
becomes spin gapfull(gapless) as has already been well 
examined by Nakamura et al. in the study of ID ex- 
tended Hubbard model. 30 We expect that this method 
is applicable to our zigzag chain model as well. 

In Fig. 4, we show the size dependence of the critical 
point t2 C /\ti\ determined by (Z a ) = 0, where we set the 
relation between the electron hopping and the exchange 
interaction as /Ki I ) 2 = Ji/Ji- For n = 1, our system 
reduces to the J1-J2 Heisenberg model and the critical 
point t2 C /\ti\ is well scaled by 1/N 2 . The extrapolated 
value |t2cAi| = 0.491 is very close to the known result of 
the Heisenberg model, J 2c / Ji = (t 2c /ti) 2 = 0.241. 31 For 
n < 1, we obtain the extrapolated value by assuming 
the same size dependence, where we set = 0.4, 

and J2/J1 = {t2/t\) 2 to correspond to the zigzag chain 
Hubbard model(t 1 -t 2 -J7 model). 

Fig. 5 shows the phase diagram on the i 2 /|ti|-n plane 
together with the result of the weak coupling theory. It 




VIM 

Fig. 5. Phase diagram on the t'2/\ti\-n plane. The solid cir- 
cles stand the extrapolated value of the critical point t2c/|il|- 
The open symbols, A, Q>> an d d indicate spin polarized 
ground state of finite size systems with S tota i/ S max = 1/4, 

Stotall ' Smax — 1/3, Stotai/ Smax — 1/2 and Stotai/ Smax — b 

respectively, where Stotai is the total spin of the ferromagnetic 
state and Smax is the maximum possible value of S to tal- The 
broken lines are the phase boundary between c2s2 and clsO ob- 
tained by the weak coupling theory. 8 The dashed lines present 
the boundary between the region of four Fermi points and that 
of two Fermi points in the nonintcracting model. 



shows that the phase boundary of the spin gap is close 
to that from the weak coupling theory for the region 
with n > 0.5 and £2/^1 1 ~ — 1- It also suggests that 
the boundary is almost independent of U in the ii-t 2 - 
U model, as has already reported for n = l. 32 ' 33 We 
note that the phase diagram suggests that the parameter 
region corresponding to Pr247 belongs to the spin gapped 
phase with Z a > 0. 

In the region 0.5 < Wl^il ~ 1-9, we find 
that the ground state of the finite system is a spin 
polarized ferromagnetic state in part as shown in 
Fig. 5, 14, 16,34 where we use the systems for n = 
4/14, 8/14, 12/14, 4/16, 8/16, 12/16, 4/18, 8/18, 12/18 and 
16/18. The region of the ferromagnetic state is discon- 
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Fig. 6. The value of v tT (Z rT ) as a function of t2/|ti| at n = 
12/18,16/18,18/18 (18-sites systems with 12,16,18-electrons), 
where we set <2 = —1 for *2 / 1 * 1 1 < —1, ti = 1 for |t2/il| < 1, 
and ti = 1 for *2 / 1 * a I > 1, respectively. For all values of t2/|ti|, 
we set Ji = 0.5t^, and J2 = 0.5t| to correspond to the t\-t2-V 
model with U = 8. The open squares stand the spin gap of Ao- 
obtained from the DMRG method for n = l. 38 



nected in respect to *2/|ti | , which is considered to be 
caused by a finite size effect. Although the phase bound- 
ary of (Z a ) = is masked by the ferromagnetic state 
in the region, we have confirmed that the sign of (Z a ) is 
positive for t 2 / |*i \ Z 1 and negative for i 2 /|£i| ~ except 
the ferromagnetic phase. 

Finally, we consider the relationship between (Z a ) and 
the value of the spin gap A CT . It is known that (Z a ) 
corresponds to the expectation value of the nonlinear 
term cos(y/8(j) a ) in the sine-Gordon model which is the 
effective Hamiltonian of ID electron system. 30 Because 
this term becomes a source producing the gap, there is 
a close relation between A CT and (Z a ). 37 We find that 
A a of the infinite system is almost proportional to the 
product of v a and (Z a ) of the 18-sites system in wide 
range of the parameter £2/1*1 1 at n — 1, where v a is 
the spin velocity corresponding to the energy scale of 
the spin part of the effective Hamiltonian. 35 In Fig. 6, 
we plot (Z a ) in the 18-sites system together with 
obtained by DMRG method for n = l. 38 

Remarkably, a phenomenological relation, A a = 
0A8v a (Z (7 }, is observed at n — 1 for all values of ^2/ |*i | - 
Assuming the same relation is satisfied even for n < 1, 
we estimate the spin gap A CT from the value of v a (Z a ) of 
the 18-sites system. To confirm this assumption, we com- 
pare our result with A CT of the t\-t 2 -U model obtained 
by the recent DMRG method. 21 As shown in Fig. 6, we 
obtain (Z a ) <~ 0.25 and v s <~ 1.5 resulting in A a <~ 0.18 
for n = 16/18 (18-sites system with 16-electrons) with 
the parameters: — 1.0, t 2 = —1.0, J\ = 0.5, and 
J 2 = 0.5. The result is in good agreement with the 
DMRG result, i.e., A a <~ 0.16 of the corresponding t\- 
t 2 -U model with |ii| = 1.0, t 2 = -1.0 and U = 8. 21 
Then, we expect that our analysis is useful to estimate 
the spin gap, even for n < 1. When we apply the above 
method to the realistic parameter region of Pr247, we 
obtain A ~ 0.0072eV for n = 16/18 with t 1 = -O.leV, 
t 2 = -0.45eV, Ji = O.OleV and J 2 = 0.15eV. In the 
case of n = 12/18, we find A CT ~ O.OlleV. These re- 



sults suggest that the order of the spin gap amounts to 
A<j <~ 100K in the realistic parameter region of Pr247 
and is larger than T c ~ 20K. 

In summary, we investigate the one-dimensional t\-t 2 - 
J\-J 2 model as an effective model for metallic CuO dou- 
ble chain of Pr247 using the numerical diagonalization 
method. The hopping parameters of electron are cho- 
sen so as to fit the <i-band from the CuO double chain 
obtained from LDA calculation and the exchange cou- 
pling energies J\ and J 2 are estimated by the known 
results of ID d-p models. In a realistic parameter re- 
gion, we show that the Luttinger liquid parameter K p is 
greater than 0.5 and the anomalous flux quantization is 
found. It suggests that the CuO double chain is respon- 
sible for the superconductivity of Pr 247 on the basis of 
the Tomonaga-Luttinger liquid theory. We also calcu- 
late the expectation value of the twist-operator Z a and 
obtain the phase boundary of the spin gap including the 
region with large value of *2/|*i|- By comparing with the 
known result of the J\-J 2 Heisenberg model, we estimate 
the value of the spin gap A a through (Z a ) and spin ve- 
locity v a , and find the spin gap becomes ~ 100K in the 
realistic parameter region of Pr247. This result is consis- 
tent with the recent NQR experiment, where (TiT) -1 is 
suppressed in the superconducting sample as compared 
with the non-superconducting sample even above T c . 
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